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ABSTRACT 

We present evidence for Milky-Way-like, low-excitation molecular gas reservoirs in near-IR selected massive 
galaxies at z ^ L5, based on IRAM Plateau de Bure Interferometer CO[3-2] and NRAO Very Large Array 
CO[1-0] line observations for two galaxies that had been previously detected in CO[2-l] emission. The CO[3-2] 
flux of BzK — 21000 at z = 1 .522 is comparable within the errors to its CO[2-l] flux, implying that the CO[3-2] 
transition is significantly sub-thermally excited. The combined CO[1-0] observations of the two sources result 
in a detection at the 3a level that is consistent with a higher CO[1-0] luminosity than that of C0[2-l]. Contrary 
to what is observed in submillimeter galaxies and QSOs, in which the CO transitions are thermally excited up to 
7 > 3, these galaxies have low-excitation molecular gas, similar to that in the Milky Way and local spirals. This 
is the first time that such conditions have been observed at high redshift. A Large Velocity Gradient analysis 
suggests that molecular clouds with density and kinetic temperature comparable to local spirals can reproduce 
our observations. The similarity in the CO excitation properties suggests that a high, MiUcy-Way-like, CO to H2 
conversion factor could be appropriate for these systems. If such low-excitation properties are representative 
of ordinary galaxies at high redshift, centimeter telescopes such as the Expanded Very Large Array and the 
longest wavelength Atacama Large Millimeter Array bands will be the best tools for studying the molecular 
gas content in these systems through the observations of CO emission lines. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: high- 
redshift — galaxies: starburst 



1. INTRODUCTION 

Daddi et al. (2008; DOS henceforth) recently reported 
the discovery of large molecular gas reservoirs inside near- 
IR selected galaxies with star formation rates SFR ^ 100 — 
150 M0 yr^' and stellar masses M* ^ IO^Mq at redshift 
z ^ 1.5. The SFR in these systems is an order of magni- 
tude lower than in submillimeter galaxies (SMGs) and QSOs 
that were also previously detected in molecular gas emis- 
sion a t high redshift (e.g.. iGreve et al.l 12005: Riecher s"et al.l 
I2OO6'). The CO[2-l] luminosities derived from IRAM Plateau 
de Bure Interferometer (PdBI) observations of two sources — 
BzK - 4171 and BzK - 21000 at z = 1 .465 and z = 1 .522 re- 
spectively — imply star-formation efficiencies (SFEs) similar 
to those of local spirals, an order of magnitude lower than 
those in submillimeter selected sources. The low SFEs sug- 
gest that these objects behave as scaled-up versions of lo- 
cal spirals, with long lasting, low-efficiency star-formation. 
Given the disk-like UV rest-frame morphology of these ob- 
jects in deep HSTh-ACS imaging, and the good agreement be- 
tween the UV and radio/mid-IR star formation rates, DOS sug- 
gested that star formation in these galaxies is taking place in 
disk-wide extended gas reservoirs. This is contrary to major 
merger-triggered nuclear starbursts as seen in local ULIRGs 
and distant SMGs, despite the fact that the SFRs and star for- 
mation surface densities in these objects are similar to those 
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of local starbursts and ULIRGs. Given that both galaxies ob- 
served in CO by DOS were detected with large CO luminosi- 
ties, and that the space density of similar sources is 10-30x 
larger than that of SMGs (Daddi et al. 2007), this suggests that 
such a quiescent gas-consumption activity could be a preva- 
lent mode of star-formation and galaxy growth in the distant 
universe. 

In order to gain insights on the physical properties of 
the molecular gas and thus on the nature of the star for- 
mation in these galaxies it is important to investigate the 
CO excitation_properties. In galaxies like the Milky Way 
dFixsen et al. 1999) and local star-forming spiral galaxies 
(e.g., Braine & Combes 1992; Crosthwaite & Turner 20()7 
Mauersberger et al. 1999; iPapadopoulos & Seaquisti 11998 



Young et al.,,1995>) the molecular gas has relatively low ex- 
citation and is rather diffusely distributed (low density and 
temperature), resulting in a high CO to H2 conversion factor. 
In contrast, dusty starburst systems like ULIRGs in the local 
Universe and SMGs and QSOs i n the distant universe have 
very dense an d warr n gas (e.g.. [Pownes & S olomonl ll998t 
Riechers et all 120061; iTacconi et all l2006l 12008; Wei B etal.l 



2005allbl) . with a CO spectral energy distribution show- 
ing thermalized^mission up to a t least the rotational tran- 
sition / = 3 (IWeifi et al.1 l2007bl) . and a CO to H2 con- 
version factor about 5 times l ower than in typical spirals 
dSolomon & Vanden Boutll2005h . Despite many CO observa- 
tions at high redshift, no evidence for Milky- Way-like CO ex- 
citation has been found so far This result may be biased by 
the pre-selection of the most IR-luminous galaxies for study, 
and by the prevalence of high frequency, high-J transition ob- 
servations. 

Here we present multi-J CO observations based on PdBI 
and Very Large Array (VLA) observations for two z 1.5 
BzK-selected galaxies previously detected in CO[2-l] with 
the aim to estimate the excitation properties of CO in these 
sources. 



2 



Dannerbauer et al. 



2. OBSERVATIONS 

2.1. Plateau de Bare Observations 

We observed BzK - 21000 at z = 1.522 in the CO[3-2] 
transition (rest-frame 345.796 GHz) with the PdBI in D- 
configuration in May and June 2008. The on-source ob- 
serving time was 9.7 hours with the full array (six antennae) 
and a beam size of 4.0" x 3.1". We tuned the receivers to 
136.970 GHz (2mm window), slightly offset to the frequency 
of the CO[3-2] line at z = 1.522 at 137.112 GHz. The to- 
tal bandwidth of our dual polarization mode observations was 
1 GHz (corresponding to ^^22001™ s^'). The correlator has 8 
independent units each covering 320 MHz (128 channels each 
with a width of 2.5 MHz). The phase center was 7.5" away 
from the position of BzK - 21000 in order to target the SMG 
GN20 during the same observations (Daddi et al. 2009). All 
the values quoted below were corrected for a primary beam at- 
tenuation (PBA) of 1 1%. We reduced and calibrated the data 
with the GILDAS softwai-e packages CLIC and MAP. The 
final data, imaged using natural weighting, have a noise of 
0.11 mJy integrated over the full 1 GHz bandwidth. In addi- 
tion to the CO[2-l] dataset already presented by D08, we also 
use CO[2-l] observations of both BzK galaxies obtained with 
IRAM PdBI in B -configuration (synthesized beam ^ 1.3") 
that have been already presented in Daddi et al. (2009). The 
new data solidly confirm the CO[2-l] detections of D08. Our 
flux calibration in both the 2mm and 3mm observations with 
IRAM PdBI is primarily based on measurements of the stan- 
dard calibrator MWC34-9 for which accurate models of the 
spectrum at these frequencies are available. The knowledge 
of the typical ranges of PdBI antennas efficiencies at the var- 
ious frequencies can also be used as a cross check of the flux 
calibration. Overall, we estimate that the typical accuracy on 
the absolute flux calibration for our observations are within a 
range of ^^10% in the 3mm band and ^22% in the 2mm band. 
To the measurement errors, we add in quadrature uncertainties 
corresponding to half those ranges for the estimate of the to- 
tal error budget on fluxes and luminosities in each observing 
transition. 

2.2. Very Large Array Observations 

We observed both sources — BzK - 2 1 000 and BzK - 4 1 7 1 
— with the VLA in D-configuration in August and Septem- 
ber 2008 (synthesized beam of ~ 1.5"). In order to mea- 
sure the CO[1-0] transition (rest-frame 1 15.271 GHz) we used 
two contiguous 50 MHz IF bands tuned at 45.6851 GHz 
and 45.7351 GHz for Bz/T - 21000 and at 46.7351 GHz 
and 46.7851 GHz for BzK -4111. The total bandwidth of 
100 MHz corresponds to about 650 km s^' in velocity space, 
which covers well most of the emission lines. The total in- 
tegration time was 21 h and 28 h on-source for BzK — 21000 
and BzK — 4171, respectively. We reduced and calibrated the 
data with standard techniques using AIPS. Given that sum- 
mer time observations with the VLA can be subject to bad 
phase stability, we employed fast switching phase calibra- 
tion using the calibrator J1302+578. We also included cycles 
where we switched to a position near to the phase calibra- 
tor, and used these data to check the phase coherence time. 
From this we measured the flux correction factor due to de- 
coherence of 24% for BzK - 2 1 000 and 1 5 % for BzK -4171. 
We correct our measurements by these factors to account for 
these small losses of signal. The final maps (100 MHz band- 
width) have rms noise levels of 0.11 mJy and 0.13 mJy for 
BzK — 21000 and BzK ~ 4111, respectively, once corrected 
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Fig. 1.— Top panel: CO[3-2] map centered on BzK-21000, aver- 
aged over 525 km s^', corresponding to the observed range of CO emis- 
sion (see also Fig. |2). Contour levels start at ±2f7 and are in steps of Iff 
(~ 0.215 mJy beam^') with positive (negative) contours shown as solid 
(dashed) lines. The white cross shows the VLA 1 .4 GHz position. Bottom 
panel: map averaged over the remaining channels. 



for de-coherence and PBA at the 10% level. For the VLA ob- 
servations, by monitoring the bootstrapped flux of the phase 
calibrator over the observing days, we estimate a gain cali- 
bration uncertainty in a range of 14%, small compared to the 
typical measurement uncertainties for our objects (see next 
section). 

3. RESULTS 

3.1. CO[3-2] in BzK -21000 

The CO[3-2] transition in BzK-21000 is securely detected 
in the IRAM PdBI data at the --6(7 level (Fig[T])- In order 
to properly compare the fluxes from the different CO tran- 
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sitions, and to derive a CO spectral line energy distribution 
(SLED), we need to perform flux measurements over the same 
velocity range in the different transitions. We use a veloc- 
ity range of 525 km s^' centered on the CO[2-l] line (the 
transition detected with the highest S/N ratio, Fig. |2]i. We 
also corrected line fluxes for possible underlying continuum 
emission. Direct measurements in the PdBI maps return no 
significant evidence for continuum emission (Fig. [T] and |2]l 
with formal measurements of 5(3mm) = 40 ± 60/iJy and 
5(2mm) = 320± 150^Jy. Given the estimate of Ljr - IO'^Lq 
in these galaxies (DOS) we would expect continuum emission 
at the level of 5(3mm) = 30/xJy and 5(2mm) = 100/xJy based 
on the SED models of Chary & Elbaz (2001) that have been 
shown to accurately reproduce the properties of BzK galaxies 
(e.g., Daddi et al. 2005; 2007). These predictions are con- 
sistent with our measurements. Given the non-detections of 
the continuum in our observations, we used the model es- 
timates as our best guesses and subtracted those values to 
obtain continuum-free line measurements. These continuum 
corrections are small, 7% and 2.5% for the 2mm and 3mm 
band, respectively. We did not correct the VLA measurements 
for continuum emission, which is expected to be entirely neg- 
ligible at 45 GHz « lO^iJy). 

Applying all the corrections and adding all sources of un- 
certainty in quadrature, our CO measurements for BzK — 
21000 ai-e /co[3-2] = 0.70 ±0.15 Jy km s"' and /co[2-i] = 
0.62 ±0.07 Jy km s^^ This new CO[2-l] measurement is 
slightly lower than that reported in DOS but it has a higher 
S/N ratio. The CO[2-l] spectrum now clearly shows a dou- 
ble peak profile, similar to BzK — 4111 in DOS, suggestive of 
rotation. 

In the following we express the line ratios in terms of 
brightness temperature ratios (or equivalently, luminosity ra- 
tios) of the different transitions, as can be derived from 
equation 3 of Solomon & van den Bout 2005): Oy-i = 
IcolJ-(J-i)]/Ico[{J-i)-{J-2)] * {{J - ^f/J^-)- In the case of 
thermalized transitions the temperature ratio is 1 by defini- 
tion, the CO luminosity is the same in the different transitions 
and the CO fluxes in Jy km s^' scale proportionally to J^. 

For BzK — 21000 we derive a solid measurement of 02 = 
0.50 ±0.12. The r32 measurement clearly shows that the 
CO[3-2] transition is sub-thermalized with respect to CO[2- 
1], a result significant at the > 4(7 level. From Fig.|2]it seems 
that CO[3-2] might be particularly weak in the reddest half of 
the spectrum, a possible sign of differential excitation. Mea- 
suring r32 for the red and blue half separately we do find a 
suggestion for different excitation (r32 = 0.2S and 0.64, re- 
spectively), but the difference is significant only at the 1.7(7 
level. 

3.2. C0[ 1-0] constraints 

Inspection of the VLA observations of the two BzK galax- 
ies shows ^ 2c7 peak signal close to the positions of both 
galaxies. Individual fluxes cannot be reliably derived from the 
data. In order to increase the S/N in the measurements we cor- 
rected each VLA map for de-coherence and PBA, aligned the 
maps at the VLA 1 .4 GHz continuum coordinates of the two 
galaxies and averaged the data from the two objects. Fig. [3] 
shows that some positive CO [1-0] signal is detected at the ex- 
pected position, significant at about 3.1(7 and corresponding 
to a peak flux density of 0.2S ± 0.09mJy. Given the insuffi- 
cient S/N of the data, it is impossible to know if the CO[l- 
0] emission is slightly extended. In order to compare this 
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Fig. 2.— The CO spectra of BzK - 21000. CO[2-l] is in tlie top panel 
(sampling of 25 km s"') and CO[3-2] is in the bottom panel (sampling of 
75 km s^'). The red horizontal line indicates the velocity range used for 
measuring the CO fluxes. The spectra are continuum subtracted as discussed 
in the text. 
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Fig. 3.— Average map of CO[1-0] VLA observations of Bz^-21000 and 
BzK — A\l\, covering 650 km s^' well centered on both Unes. Contours start 
at ±2a. The map has a noise of 90/JJy. The cross coiTesponds to the VLA 
1.4GHz continuum positions. 
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Fig. 4.— Left panel: The CO SLED of Bz^T- 21000 at z = 1.522 obtained from PdBI a nd VLA observations. The black solid, dotted and dashed lines show 
our LVG models (model 1, 2, 3, respectively). The Milky Way SLED from lFixsen et alj (triangles. 1199 9) is shown for comparison. The constant brightness 
temperature relation (red dashed line) is normahzed to our CO[2-l] measurement. Right panel: The same of left panel but shown in luminosity units, normalized 
to the C0[2-l] transition. Juppe,- > 2 transitions can severely underestimate the CO luminosity and thus the molecular gas mass. 



CO[1-0] measurement with the CO[2-l] results we have used 
the B -configuration CO[2-l] observation with IRAM PdBI, 
which have comparable angular resolution to the VLA CO[l- 
0] data, and extracted the signal over velocity ranges match- 
ing the VLA observations on each galaxy (covering very well 
both lines). This results on an average peak CO[2-l] signal 
in the combined CO[2-l] map of 0.67 ± 0.09 mJy. We derive 
r2i = 0.60^q |7, averaged over the two BzK galaxies. Taken at 
face value, this implies that we are seeing excess peak emis- 
sion in CO[1-0] relative to what is expected from CO[2-l] 
although the evidence for an excess is very mild given the low 
S/N in the CO[1-0] data. Higher S/N observations would be 
needed to establish if the CO[2-l] transition is sub- thermally 
excited with respect to CO[1-0]. 

4. DISCUSSION 

For BzK -21000, the CO[3-2] and [2-1] fluxes are com- 
parable within the errors. The CO spectral line energy dis- 
tribution of BzK — 2100 is consistent w ith the Milky- Way 
disk measurements by Fixsen et alJ (1199 9) (Fig.lH), and simi- 
lar to what is observed in local spirals (e.g ., Braine & Combea 
1992'; 'Crosthwa ite & TurnedllOOTi: iMauersberger et al., ,19991 
Young et al. 1991). The CO SLED of BzK - 21000 is in 
stark contrast to what has been generally seen so far in high- 
redshift sources like QSOs and SMGs, for which gas has 
been f ound to be thermalized up to at least the ■/ = 3 tran- 
sition (iRiechers et al.ll2006l: IWeiB et alJl2005bl l2007allbh . A 
possible exception is the extrem ely red galaxy HR 1 at 
z = 1 .4 4 according to the results of' Andreani et al. (2000) and 
iPapadopo ulos & Ivison (20i)2) who find comparable fluxes of 
CO[5-4] and CO[2-l]. However we note that the excitation in 
our sources is still lower than that in HRIO. 

We now discuss the implication of our excitation mea- 
surements on the physical properties of the gas, includ- 
ing the kinetic temperature {Tki„) and density (n) of the 
H2 gas in the molecular clouds , by using Large "Veloc- 
ity Gradient ( LVG ) modeling (IGoldreich & KwanI 119741: 
Scoville & Solomo5[i97l. We use the colUsion rates from 
Flowed (12001) with a fixed CO abundance to velocity gradient 
ratio o f [CO]/(dv/dr)= 1 x lO^^ pc (km s'^)-^ (WeiB et alj 
l2005ah . an H2 ortho:para ratio of 3:1 and the cosmic back- 



ground temperature at z = 1.522 of 6.89 K ^. Due to the 
limited amount of information available, we will not be 
able to resolve the ambiguity between all the free parame- 
ters. E.g., choosing a different value of [CO]/(dv/dr) would 
also affect the SLED, and would require re-adj usting T^m 
and/or n to reproduce observed ratios (see, e.g., IWeiB et alj 
2005a, for a more detailed discussion of the parameter space 
degeneracies). We present here three representative mod- 
els, spanning the range of parameter degeneracies in Tj;,,,, 
and n (see Fig. Hjl, that reproduce the observed CO SED 
of BzA' — 21000. Model 1 has a kinetic gas temperature 
Tkin = 25 K and molecular cloud density of n = 1300 cm^^, 
very similar to the molecular gas excitation conditions of 
the Milky-Way. The filling factor of clouds, assuming a 
disk radius ro ^ 5 kpc from our ACS observations (DOS), 
is 2% for this model. Model 2 can still fit the data assum- 
ing 7i,„ = 90 K and n — 600 cm^^ with a similar filling fac- 
tor to model 1 . Model 3 results in a higher filling factor of 
8% with a low = 10 K but a higher n = 2500 cm^^ 
Typical gas/dus t temperatures in nearby spirals/starbursts are 
10-30 K (e.g., i Braine & Combeslll992l: IMauersberg er et alJ 
Il999l) . sometime s up to 40-50 K in the cent ral regions of nu- 
clear starbursts ( iDownes & Solomonlll998l) . Normal galac- 
tic giant molecul ar clouds have typical average den sities of 
10^ - 10^ cm-3 dEvanslll999l: IScoville et alJri974h . These 
values are broadly similar to what is allowed from LVG mod- 
eling of our observations. Clearly, we can rule out solutions 
where the gas is at high density and/or temperature, similar to 
what is generally found inside local ULIRGs, distant SMGs 
and QSOs. 

Despite the degeneracies in the inferred physical param- 
eters, the LVG models presented above consistently predict 
r2i ^ 0.85, higher than our average estimate for the BzK 
galaxies but consistent within the uncertainties. Robustly con- 
firming an excess CO[1-0] emission over what is predicted by 
the LVG models shown in Fig. |4] would imply the discovery 

' Even for the lowest CO[1-0] transition the cosmic background radiation 
at z= 1.5 will have only a minor effect on the depopulation of the state of 
our two sources. The reason is that the intrinsic temperature in the molecu- 
lar clouds is expected to be higher, e.g. the intrinsic temperature within the 
Milky- Way clouds is at least greater than 10 K and most likely lies around 
~20 K. 
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of very cold molecular gas unseen even in the CO[2-l] transi- 
tion. For the moment, given our limited S/N at the measured 
r2i, we assume that r2i ^ 0.85 and r32 ~ 0.50 are the typical 
values for our sources. This implies that the CO luminosities 
derived for our galaxies in DOS from CO[2-l] do not need 
substantial revisions, exception made for the small 16% in- 
crease required when converting CO[2-l] to CO[1-0], as sug- 
gested by our LVG modeling. However, our results strongly 
indicate that when using transitions higher than CO[2-l] only 
a fraction of the CO[1-0] luminosity is recovered (Fig. IH 
right), implying that such measurements are not reliable for 
accurately estimating the total gas mass and star formation ef- 
ficiencies in BzK-like galaxies. This evidence could explain 
the discrepancy between the observations of Hatsukade et al. 
(2009) and Tacconi et al. (2008) who failed to detect CO[3- 
2] in 4 optically and UV selected galaxies at 1.4 < z < 2.5, 
while so far we have a high detection rate in CO[2-l] at red- 
shift z 1.5. A large correction factor of 2.4 is required, in 
fact, at CO[3-2] in order to recover the total molecular gas 
mass, based on our results from BzK — 21000. 

Given that these two high redshift galaxies have star- 
formation efficiencies and CO excitation properties similar to 
those in local spirals and in the Milky Way, it is plausible that 
they may also have a Milky-Way-like molecular conversion 
factor Oco as well. This would imply molecular gas reservoirs 
of ~ 10" Mq and gas mass fractions fg^s ^ 0.6. To summa- 
rize, our sources provide clear examples of the long-sought 
after population of low-excitation CO emitters at high-z. Ob- 
servations of a larger sample of near-IR selected galaxies are 
required in order to confirm our results in a statistical way. In 
particular, high S/N observations of CO[1-0] and high-J CO 
are required to search for additional components of colder and 
warmer gas, respectively. Much higher S/N data will also be 
required to search for excitation variations inside the disks, as 
tentatively suggested by our data and as could be expected in 
the case of very clumpy gas distribution (e.g., Bournaud et al. 
2008). 

Given that these are the first excitation measurements in- 



side normal, near-IR selected disk-like galaxies in the dis- 
tant universe, it is plausible to speculate that such low- 
excitation properties might be typical in distant galaxies. 
ALMA will be a powerful device to study CO[2-l] in z < 3 
galaxies, but at much higher redshift, the higher order tran- 
sitions may not be excited i n the average galaxy (see also 
iPapadopoulos & Ivisonl20()2h . In those cases, centimeter tele- 
scopes, such as the EVLA or the Square Kilometer Array 
(SKA), will become the primary tools for study of molecular 
gas in the earliest, normal galaxies. In distant galaxies with 
massive amounts of cold, diffuse, low-excitation gas we will 
be able to detect carbon lines with ALMA (e.g. [CII]A 158 jUm 
redshifted to > 400/im in z > 1.5 galaxies). This synergis- 
tic combination of observations will offer a powerful tool for 
interpreting the gaseous content of ordinary galaxies at high 
redshift. 
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